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Short, hard, 7-ray bursts (GRBs) are believed to originate from the coalescence of two neutron 
stars (NSs) or a NS and a black hole (BH). If this scenario is correct, short GRBs will be accompanied 
by the emission of strong gravitational waves (GWs), detectable by GW observatories such as LIGO, 
Virgo, KAGRA, and LIGO-India. As compared with blind, all-sky, all-time GW searches, externally 
triggered searches for GW counterparts to short GRBs have the advantages of both significantly 
reduced detection threshold due to known time and sky location and preferentially higher GW 
amplitude because of face-on orientation. Based on the distribution of signal-to-noise ratios (SNRs) 
in candidate compact binary coalescence (CBC) events in the most recent joint LIGO- Virgo data, 
our analytic estimates, and our Monte-Carlo simulations, we find an effective sensitive volume for 
GRB-triggered searches that is «2 times greater than for an all-sky, all-time search. For NS-NS 
systems, a jet angle 6j = 20°, a 7-ray satellite field of view of 10% of the sky, and priors with 
generally precessing spin, this doubles the number of NS-NS — short GRB and NS-NS — short GRB 
associations, to ~3-4% of all detections of NS-BH and of NS-BH. We also investigate the power 
of tests for statistical excesses in lists of sub-threshold events, and show that these are unlikely to 
reveal a sub-threshold population until finding GW associations to short GRBs is already routine. 
Finally, in Appendix A we provide useful formulae for calculating the prior distribution of GW 
amplitudes from a CBC, for a given GW detector network and given sky location. 



I. INTRODUCTION 

We currently sit between the first and second gener- 
ations of kilometer-scale, ground-based interferometric 
GW detectors. The first direct detection of GWs will 
very like occur before the end of the decade. 

The first detected signals will probably be from merg- 
ers of neutron star-neutron star (NS-NS), neutron 
star-black hole (NS-BH), and black hole-black hole 
(BH-BH) binaries, collectively referred to as compact bi- 
nary coalescences, or CBCs. NS-NS and NS-BH coales- 
cences are also likely progenitors for most short GRBs [ , 
] . It is therefore natural to use detections of short GRBs 
to trigger searches for GW signatures of CBCs that oc- 
cur at the same instant (to within a few seconds) and 
the same sky location (within the error bars). Such 
GRB-triggered searches for GWs from CBCs are already 
being carried out [3-5]. 

In this paper, we address several questions regarding 
GRB-triggered CBC searches. We begin by reviewing 
recent results from such searches in Section I B and con- 
sider what GW and GRB detectors will be available in 
the Advanced detector era in Section II. In Section III, we 
discuss jet collimation and the reduced detection thresh- 
old from knowing the GRB time, which are the biggest 
factors that play into detection prospects. Our primary 
results appear in Section IV, in which we estimate the 
rate of coincident detections both analytically and via 



Monte-Carlo simulations. A closely related data analyis 
activity has been the search for a statistical excess of 
high-SNR, but sub-threshold, CBC candidate events co- 
incident with short GRBs. In Section V we predict the 
science yield from searching for such statistical excesses, 
and demonstrate that the extra information will typically 
be negligible. Finally, in Appendix A we derive a num- 
ber of useful analytic formulae for describing a detector 
network's sensitiviy to a CBC at a given sky location (in 
the case of interest here, the location of the GRB). 

We note that rates of short GRB-GW coincident de- 
tections were also estimated in recent papers by Chen 
and Holz [6] and another by Kelley et al. [7] which ap- 
peared when this paper was almost finished. Our method 
is similar to both of theirs. Our conclusions are qualita- 
tively similar to those of [6], but are qualitatively differ- 
ent from [7] due to different assumptions and approxi- 
mations: [6, 7] both assume a Gaussian distribution of 
outliers, while we base our calculations on the distribu- 
tion of outliers observed in Laser Interferometer Gravita- 
tional-Wave Observatory (LIGO) 's sixth and Virgo's sec- 
ond and third science runs (S6/VSR2,3). [7] also adopt 
a fiducial value for the GRB jet angle that is a factor 
~ 2 smaller than ours. (Both values are plausible, given 
the uncertainties.) Because of our different assumptions 
about the statistics of the GW search and the GRB jet 
angle, we derive a GRB-GW detection rate that is ~ 20 
times higher than theirs, which lifts it from ~ 1/10 yr to 
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~ 2/yr (i.e., from almost negligible to interesting). 



A. Science motivation 

The detection of GWs from CBCs will have several 
scientific implications. The masses of the two compact 
objects will be determined quite accurately [8, 9]. With 
enough SNR, the spins of these objects can also be con- 
strained [10]. These measurements, and the overall rates, 
will provide information on stellar evolution [11]. The 
details of the late inspiral and post-merger gravitational 
waveform will also inform the high-density NS equation 
of state [12, 13]. Details of the merger will also permit 
tests of general relativity in the strong-field regime [14], 
tests of local Lorentz invariance [15], and constraints on 
the graviton mass [16, 17). 

A coincident short GRB-GW detection would prove 
that at least some GRBs are indeed produced by merger 
events. Furthermore, it should be possible to determine 
the redshift of the short GRB's host galaxy, while the 
GWs accurately encode the distance to the binary. It 
has been shown that 10 short GRB with redshift mea- 
surements could constrain Hq to within 2% assuming a 
GRB jet angle of 20° [18]. 



B. Recent results 

To date, two types of searches for CBC signals associ- 
ated with short GRBs have been executed: single-event 
targeted analyses for short GRBs associated with very 
nearby galaxies [GRB 070201 and GRB 051103; 3, 19], 
and analyses covering all short GRBs during LIGO and 
Virgo data-taking epochs [5, 20]. None of these analyses 
found a significant GW candidate, but the results were 
used to establish lower limits on the distances, assuming 
CBC progenitors. 

GRB 070201 and GRB 051103 had localizations that 
significantly overlapped with the galaxies M31 [An- 
dromeda, 770 kpc away; 21, 22] and M81 [3.6 Mpc away; 
'•'!], respectively. However searches in contemporaneous 
GW data were able to exclude CBCs as the source of 
those bursts [3, 19]. 

The 22 short GRBs that occurred during the 
LIGO-Virgo data-taking run LIGO's fifth and Virgo's 
first science runs (S5/VSR1) were followed up with CBC 
searches in the GW data, using analysis methods similar 
to those for GRB 070201 and GRB 051103 [5]. The 26 
short GRBs that occurred during S6/VSR2,3 runs were 
also followed up [20], using an improved, coherent analy- 
sis strategy [24]. No coincidences were found, and lower 
distance limits on putative CBC counterparts were estab- 
lished. In the more sensitive, S6/VSR2,3 search, the me- 
dian 90% confidence lower limits for NS-NS and NS-BH 
binaries were 16 and 28 Mpc, respectively. Both analy- 
ses included a test for a sub-threshold population excess. 
For S5/VSR1, a Mann- Whitney U test [25] was used, 



while a binomial test was used for S6/VSR2,3. In both 
cases, the sub-threshold populations were found to be 
consistent with the background. 

II. DETECTOR NETWORK ROADMAP 

To make sensible predictions of the outcome of future 
GRB-triggered GW searches, one needs to know what 
GRB and GW detectors might be operating in the next 
decade. 



A. GW Detector network roadmap 

The US LIGO [ ] has recently completed a one-year 
data-taking period between July 2009 and October 2010, 
in coincidence with the French-Italian Virgo detector 
[27]. 

LIGO is currently upgrading to its Advanced detector 
configurations [28, 29], with the goal of increasing the 
sensitivity gradually to a factor of ~10 compared to the 
initial configuration and extending seismic-limited sensi- 
tivity to lower frequencies. More recently, the US Na- 
tional Science Board has authorized one LIGO detector 
to be moved to India in order to vastly expand the world- 
wide detector network's sky localization capabilities. The 
earlier attempts to do the same in Australia have been 
formally abandoned. 

Virgo is upgrading to the Advanced Virgo configura- 
tion, similar to Advanced LIGO in optical layout and 
sensitivity. The start of the data-taking period with the 
Advanced detectors is foreseen in ^2015. 

GEO600 is a British-German detector with 600 m arms 
and Advanced optical configurations [30]. GEO600's 
next-generation configuration will be GEO-HF with a 
focus on high-frequency sensitivity [31]. GEO-HF 's 
high-frequency sensitivity will be the best in the world, 
which will be useful in parameter estimation [ ]. How- 
ever, low-frequency sensitivity is more important for 
CBC detection and GEO600's relatively short arms put 
it at a significant disadvantage. It will likely not be used 
for detection searches. 

Construction has begun on KAGRA (KAmioka GRAv- 
itational-wave observatory; formerly LCGT) in Japan 
[32], which should debut in 2018 with its design sensi- 
tivity. It has 3 km- long arms constructed underground 
and uses cryogenically cooled sapphire mirrors for test 
masses. The final detector is expected to detect a NS-NS 
system at a distance of 240 Mpc with SNR=10 [33]. 

B. GRB Detector network roadmap 

During the last S6/VSR2,3 science run, the triggers 
came mostly from the Swift [34] and Fermi missions [35], 
with a few from the Interplanetary Network (IPN) [36]. 
It should be noted here, that IPN also detected most of 
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the Swift /Fermi triggers, but with a much poorer sky 
location because only triangulation methods can be used 
by IPN. 

IPN has unfortunately lost its primary funding, but 
is nonetheless expected to operate during 2015-2020, if 
perhaps with a smaller number of satellites, still detect- 
ing GRBs at these times but with a lower rate [37]. 
Swift might operate for another 5 years or even longer, 
since it has no expendables and the spacecraft is in good 
shape [38], but the operation depends on NASA funding. 
The Fermi instruments might also operate until ~ 2018. 
The 7-ray Burst Monitor (GBM) instrument on Fermi 
achieves instantaneous sky coverage of about 70% or 
8.8 sr [30% of the sky being occulted by the Earth at 
the altitude of Fermi's orbit; 39], but GRBs detected by 
the GBM alone are very poorly localized. 

Lobster is a proposed NASA mission similar to Swift in 
strategy, with a wide-field X-ray imager (WFI), narrow- 
field followup IR telescope (IRT), and slewing apparatus 
to point the latter [40]. WFI is more sensitive than the 
Swift Burst Alert Telescope (BAT), but has a smaller 
field of view (FOV) at 0.5 sr. 

The French-Chinese Space-based multi-band astro- 
nomical Variable Objects Monitor (SVOM) mission is 
targeted at a broader scientific target, including answer- 
ing questions related to GRBs, cosmology, and funda- 
mental physics [ : ]. Its main instrument, ECLAIR, is a 
coded-aperture telescope aimed at a broad energy range 
of 4-250 keV, with a FOV comparable to that of Swift. 
The effective detection area is also close to that of Swift, 
resulting in an expected number of 70-90 GRBs per year, 
of which 20-25% could be short. As the ECLAIR tele- 
scopes are more sensitive to a lower energy range com- 
pared to BAT, the extended emissions and afterglow of a 
GRB can be observed more deeply, resulting in improved 
GRB locations, and hence in a larger number (—50%) of 
GRBs with redshift measurements [42]. The anticipated 
launch date is -2015-2020. 

The South Korean-led Ultra-fast Flash Observatory 
Pathfinder (UFFO-P) mission intends to catch the rise 
of GRBs [43]. It has been constructed and is anticipated 
to launch in June 2013. It carries a coded-aperture burst 
alert telescope similar to Swift's BAT, sensitive from 15- 
200 keV and with a FOV of -2 sr. UFFO-P's headlin- 
ing feature is that it can repoint in response to a trig- 
ger in — 1 s using its slewing mirror telescope (SMT) , 
which is a substantial improvement in response time over 
Swift's — 1 minute to slew the whole spacecraft. Though 
UFFO-P has a small collecting area and only a small 
optical telescope for followup, the pathfinder already has 
some discovery potential. The conceived UFFO-100 mis- 
sion will increase collecting area, replace the SMT with a 
still faster MEMS micromirror array to redirect its opti- 
cal path and add an NIR camera, with the goal of gather- 
ing a statistically significant population of rising GRBs. 

In conclusion, several missions are expected to oper- 
ate during the Advanced detector area, which are ei- 
ther already operating (like Swift, Fermi and IPN3), are 



in development (like SVOM and UFFO-P), or planned 
(like Lobster and UFFO-100). However, given the un- 
certainty in how many of these missions will ultimately 
fly, throughout this paper we will assume that during the 
Advanced GW detector era the effective coverage of the 
combined GRB detector network will be approximately 
that of Swift's BAT alone, 1.4 sr [ ] or about one-tenth 
of the sky. 



III. DETECTION PROSPECTS 

A. Collimation 

GRBs show strong evidence for collimated, relati vis- 
tic outflow along a jet. Assuming that the jet is roughly 
conical, its size is described by the jet angle 0- 3 , from the 
center to its outer edge. We define /& to be the fraction 
of the sky into which 7 rays are launched. For a CBC 
that emits a single jet, this is / b = (1 — cos6>j)/2. If the 
CBC emits two jets, presumably in opposite directions, 
then this fraction doubles to /& = (1 — cosflj). Collima- 
tion reduces the number of CBCs that are observable as 
short GRB events, since the observer only sees the GRBs 
for which the Earth lies within the jet. 

The Lorentz factor Y of the beam decreases as it sweeps 
up external material, and at the point where Y reaches 
— 1/0 j, the flux decay abruptly steepens due to special 
relativistic effects. The beaming half-angle can be deter- 
mined for a given GRB by the time of this "jet break." 
However the rapid decay of the late-time lightcurves of 
short GRBs makes the estimation of 6j difficult. Grupe 
et al. [ ] places a lower limit of 6*j > 25° while Bur- 
rows et al. [46] infers the value of 9j to be in the range 
4-8°. Goldstein et al. [47] suggest a value in the range 
between 40° and 90°. Fong et al. [48] find 3-8° in a re- 
cent short GRB. Simulations of NS-BH mergers indicate 
a range of 0- 3 ~ 30-50° for binaries of moderate spin while 
finding 5-10° for near-extremal spinning systems [49]. In 
this paper we will adopt 0j — 20° as a fiducial value when 
quoting results; however our analytic estimates and Ta- 
bles allow the reader to trivially convert the results to 
other values of y 

Of course, the GRB beaming angle is presumably not 
a universal number — it has some distribution function. 
While our simulations assumed that all GRBs had the 
same value for /;,, our results on rates in this paper are 
approximately generalizable by simply replacing by 
its average value (/{,). We emphasize that (/&) refers to 
the average over all short GRBs, not just the detected 
ones, since the detected population depends on selection 
effects. (E.g., the detected population is biased towards 
GRBs with larger values of at fixed flux.) 
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B. Reduced search space 

The search for a GW CBC signal triggered on an 
electromagnetic (EM) counterpart has sensitivity advan- 
tages over an all-sky search. Semi-analytic calculations 
and numerical simulations predict that the majority of 
the NS matter is accreted within milliseconds to seconds, 
so it has become common practice to search for GW trig- 
gers within a [—5, 1] s 'on-source window' surrounding a 
GRB to account for up to a 5-second GRB-GW delay and 
up to 1 s of uncertainty in the GRB time-of-arrival (TOA) 
[see sections 2.2 and 5.1 of 20, for references]. There are 
further possible reductions in the search volume due to 
the known sky location and even by restricting to the 
space of CBC parameters that allow for tidal disruption 
outside the innermost stable circular orbit, but in this pa- 
per we will neglect them, as they have much less impact 
than the reduction in observation time. 

In this section, we are interested in estimating the re- 
duction in SNR threshold in going from an all-sky search 
to a GRB-triggered search while keeping constant the 
false-alarm probability (FAP) at the detection threshold. 
The first detection is likely to be held to a high stan- 
dard of FAP < 10~ 6 , but once detections are routine, 
the threshold should be lessened considerably. Through- 
out this paper, we assume that FAP < 10 -4 is required. 

For 20 short GRBs per year of livetime (T a n = 1 yr) , 
the observation time for GWs is T gr b = 120 s. Assum- 
ing that the searches are comparably effective at back- 
ground rejection, the false-alarm rate (FAR) at a given 
SNR should be the same, but the FAR at a given SNR, 
FAP(p) = TFAR(p), is reduced by a factor 

FAP grb (p) T grb ^ 
FAP aU (p) T all ~ 

We estimated the SNR threshold for a GRB-triggered 
search using the background SNR distribution from the 
S6/VSR2,3 all-sky search [50]. The all-sky search was a 
coincident search, in which matched filtering and thresh- 
olding were performed on individual detectors and spuri- 
ous triggers were vetoed by demanding consistent TOAs 
in multiple detectors. For a targeted, GRB-triggered 
search, a coherent search is possible in which threshold- 
ing is done on the suitably time-delayed and summed 
SNR for the whole network [ ]. False-positive rejec- 
tion is aided by tests on null streams, which are special 
linear combinations of the detectors that are insensitive 
to GWs. Coherent searches are less feasible for all-sky 
searches because a different time-delay is required for 
each sky location. Since in Gaussian noise a coincident 
search with a two-detector network has the same statis- 
tics as a coherent search with any network of (more than 
one) detectors [24], we can extrapolate the threshold for a 
targeted, coherent search from the statistics of an all-sky, 
coincident search. 

For the S6/VSR2,3 all-sky search, Figure 3 of [50] and 
the data available online at [ ] give the FAR as a func- 
tion of p c , a x 2 -weighted quadrature sum of the SNR in 



all of the detectors that has been found to be a useful 
detection statistic. We assume that this p c is equiva- 
lent to the network coherent SNR p, as argued above. 
From this Figure, we find that during S6/VSR2,3, FAR = 
10~ 4 yr" 1 when p c = 11.3. When FAR = KT 4 /120 s, 
p c = 9. From this, we take pth.aii = 11-3 as the SNR 
threshold for an all-sky search and pth.grb = 9 as the SNR 
threshold for a GRB-triggered search. A triggered search 
can see 11.3/9 times farther than the all-sky search. On 
the other hand, as we will see, the increased range is 
somewhat offset by jet collimation and the limited FOVs 
of high-energy satellites. 



IV. RESULTS 

A. Analytic estimates 

Here we provide some simple estimates of the CBC de- 
tection rate we expect for triggered searches, compared to 
the CBC detection rate for untriggered, all-sky searches. 
The method is the same as in Kochanek and Piran [52], 
but our inputs and conclusions are different. We will as- 
sume that short GRBs come only from CBCs, and denote 
by /qc the fraction of CBCs that produce short GRBs. 
For the CBCs that produce bursts, let /b be the aver- 
age solid angle into which 7 rays are launched. Again, 
the average is over the whole GRB population, not just 
population of detected GRBs (which is biased towards 
GRBs with wider beams at fixed flux). For any CBC 
that produces two beams (presumably in opposite direc- 
tions) , we consider /b for that burst to be the sum of the 
solid angles for the two beams. Finally, we define S to be 
the average fraction of the sky that is "covered" by the 
then-existing GRB detector network. For the Advanced 
GW detectors, most CBCs will be detected at distances 
< 200 Mpc, while most detected short GRBs are much 
farther away (D ~ 1 Gpc) , so the GRB accompanying an 
observed CBC should be detectable as long as the Earth 
lies within the beam, and the source is within the tele- 
scope's FOV. Then the fraction F of CBCs for which we 
detect a GRB is 

F = f GC Sf h (1) 

We shall adopt F = 6 x 10~ 3 as a reasonable fiducial 
number. This would correspond, e.g., to S — 0.1, and all 
CBCs emitting two jets of 6j = 20° each. 

As explained above, by looking for NS binaries that are 
roughly coincident in time with observed short GRBs, we 
increase the sensitivity of the search at fixed FAR. Let 
Pth.aii/pth,grb be the ratio of the GW detection thresholds 
for blind and GRB-triggered searches, respectively. As 
explained in Section IIIB, we adopt MS M fiducial value 
Pth,all/Pth,grb ~ 11.3/9.0 « 1.25 . 

The final factor that we need arises because the merg- 
ers that we see in 7 rays are presumably the ones for 
which we are viewing the binary nearly face on (because 
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the 7 rays are presumably beamed perpendicular to the 
orbital plane.) The quadrupolar pattern of the emitted 
GWs is also strongest along the direction perpendicu- 
lar to the orbital plane; the amplitude of h on the axis 
is fa 1.51 times stronger than the isotropic detection- 
averaged value. (It is well know that that rms enhance- 
ment is y/5/2 — 1.58, but what is relevant for determin- 
ing rates is the cube root of mean-cubed enhancement, 
which is 1.51.) Although 1.51 is the enhancement factor 
of an optimally oriented binary, we show in Appendix A 
that as long as the beam half-angle is < 25°, then using 
1.51 is at most a ~ 5% overestimate, which is acceptable 
for our purposes. 

For Advanced LIGO, we will be detecting CBCs in 
the range ~ 50Mpc to ~ 500 Mpc. At these distances, 
we can safely approximate spatial geometry as Euclidean 
and the density of mergers as uniform. Let 7v a ii be the 
CBC detection rate of blind, all-sky Advanced GW de- 
tector searches, and let Af gY b be the rate for triggered 
searches. What is the ratio A/" g rb/-A/"aii? The volume in 
which GRB-triggered CBCs are detectable is larger than 
for the blind case by the factor (1.51p t h,aii/pth,grb) 3 j but 
recall that only a fraction F of the CBCs emit detectable 
7 rays. Thus 



all 



PthA 



F ~ 0.041. 



(2) 



In other words, of the first ~ 25 detections, we would ex- 
pect only one to come from the GRB-triggered pipeline; 
it is correspondingly unlikely that the very first GW de- 
tection of a CBC will result from a GRB trigger. This 
is the same conclusion reached in Chen and Holz [(>] and 
Kelley et al. [7]. 

Of the CBC detections associated with short GRBs, 
a fraction (pth,aii/Pth,grb) 3 ~ 0.5 will be strong enough 
to be detectable even without the GRB trigger. So the 
increase in the rate of GW detections, thanks to GRB 
triggers, will be ~ 2%. This is a small increase in CBC 
detections, but, given the extra information to be gleaned 
from the EM counterpart, a non-negligible one. Note 
that this 2% increase is about 2 times higher than the 
fiducial increase reported in Kelley et al. [7]. The differ- 
ence comes mostly from two factors: i) a factor ~ 4-5 
from our use of non-Gaussian statistics, and ii) a fac- 
tor ~ 6-7 from our wider fiducial beam angle of 20° vs. 
their 11°, combined with our assumption of two oppo- 
site-pointing jets per CBC. 



B. Effective sensitive volumes 

For a given detector network and mass pair, we would 
like to compute the relative detection capability of a tar- 
geted short GRB search compared to an all-sky search. 
Following Finn and Chernoff [ ], but generalizing from a 
single detector to a network, we define an effective sen- 
sitive volume I4ens, which can be multiplied by a con- 



stant rate density 1Z to obtain the total detection rate 
N = KV scns : 

V scns = Q &(p - Pth )dV^ 




Ah 



r 2 dr 



(3) 



where is the Heaviside step function, (a, 5) are the 
source's right ascension and declination, and D t h is the 
distance at which the network registers SNR above its 
threshold p t h for a given source type. depends on 

the full specification of a waveform plus the detector noise 
spectrum and the various projection angles. Eq. 3 is quite 
convenient for Monte-Carlo integration. See Appendix A 
for network SNR expressions and analytical evaluation of 
equation 3 in the non-spinning case. 

To take into account the imperfect duty cycle of the 
detectors, along with our requirement that at least two 
detectors be "on" to claim a detection, it is useful to de- 
fine a double detection volume (2DV), in close analogy to 
the triple detection rate (3DR) of Schutz [ ] . (Our for- 
mulation generalizes the treatment in Schutz [53] since it 
applies to spinning binaries and removes the assumption 
of equal detector responses.) For a set of detectors Af, 



2DV = 



E 



/l*l(l-/)l"H*IV. 



X 

sens i 



(4) 



where V>2{M) is the set of subsets of AT with size 2 or 
greater and \X\ denotes the size of set X. f is the duty 
cycle, which we take to be the same for all detectors. 



C. Results from Simulations 

We evaluate the 2DV (equation 4) via Monte-Carlo 
integration, randomly drawing the sky location and the 
polarization angles, for both NS-NS (1.4-1.4 M ) and 
NS-BH (1.4-10 M e ) pairs. We use the Taylor T4 wave- 
form [ ], including general spin precession, to determine 
expected signal amplitudes and thus D^. We match 
the distributions of spin amplitude and spin orientation 
used in the recent GRB 051103 analysis [19] with NS di- 
mensionless spin drawn uniformly from [0,0.4], BH spin 
drawn uniformly from [0,0.98], and uniform spin ori- 
entations subject to a cut on the tilt angle (angle be- 
tween BH spin axis and NS orbital axis) to be <60°. 
We also adopt their rigid rotation procedure so that it 
is the total angular momentum that is uniformly dis- 
tributed within the cone of 9j rather than the orbital 
angular momentum. Using different values for 9j, we 
computed the relative improvement in effective sensi- 
tive volume 2DV of a triggered search (restricted incli- 
nation angle, coherent threshold pth,grb = 9.0) compared 
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to an all-sky search (unrestricted inclination, coincident 
threshold pth.ail = 11.3). 

The 2DVs have been defined such that the detection 
rate for the all-sky search is A4n = TZVqqo 3 , where 1Z is 
the merger rate density in the local universe and V^' 3 the 
2DV assuming a coherent SNR threshold of 11.3 and a 
source population whose inclination is unrestricted (#j = 
90°). Similarly, the number of GRBs associated with a 
GW detection is given by A/'grb = S fh H Vg. . The rate 
for GRBs independently detected by both searches (i.e., 
with GW network SNR > 11.3) is Af hoth = S f h K V 9 f 3 . 
Our computed sensitive volumes are shown in Table I. 
Each value of V g 9 represents 10 6 simulations, while the 

Vg 1 ' 3 values are simply rescaled from Vg.. 

Under the assumption that all CBCs produce 
short GRBs and vice versa, we translate sensitive vol- 
umes to detection rates in Table II and Figure 1. We 
take values of 0.0198 Lio/Mpc 3 and "realistic" coales- 
cence rates of 6 x 10~ 5 L7/ yr -1 for NS-NS sources and 
2 x 10~ 6 Liq yr- 1 for NS-BH, as used in Abadie et al. 
[55]. The coalescence rates are uncertain by two orders of 
magnitude, but our relative volumes have much smaller 
errors, so we show two or more significant figures not to 
reflect uncertainty but rather to allow for relative com- 
parisons. We have assumed an effective FOV of S = 0.1, 
independent duty cycle factors of / = 0.8, and a net- 
work of LIGO Hanford, LIGO-Livingston, Virgo, KA- 
GRA [ ], and LIGO-India [HLVKI; orientation guess 
via 53]. For all GW detectors, we used publicly avail- 
able projections of noise spectra at full design sensitivi- 
ties [57-59]. However, we started SNR integration from 
a lower frequency limit of 40 Hz in order to fit waveforms 
within the memory limits of the bulk of computers avail- 
able to us, enabling much greater parallelism, at some 
cost in realism. These assumptions lead a total annual 
detection rate of 61 NS-NS signals and 18 NS-BH signals 
for the all-sky search. 

Assuming an effective FOV of S = 0.1 for the GRB 
detectors, and that every CBC produces two opposite 
jets, each with a fiducial jet angle of 0j = 20°, we find 
that 1.8% of CBCs detected by the all-sky search would 
be coincident with observed short GRBs. Adding the 
triggered search increases the number of coincidences to 
3.8%. Clearly, our Monte-Carlo results are in very close 
agreement with the analytic estimates in Sec. IV A. 
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V. SUB-THRESHOLD STATISTICAL EXCESS 
TESTS 

In a calendar year, there will be some number of 
searches for GWs associated with short GRBs. Some of 
these may result in individual detections, meaning that 
there are events that survive all vetoes and have SNR 
above some theshold p\ for confident single-source detec- 
tion. However one can also consider some lower threshold 
p c , and look for a statistical excess of events with SNR 



FIG. 1. Fractional number of additional CBC detections con- 
tributed by an externally triggered search (solid) and fraction 
of all-sky CBC detections predicted to be in coincidence with 
short GRB observations without a special triggered search 
(dashed). NS-NS and NS-BH ratios are equal within statis- 
tical uncertainty, so only NS-NS results are plotted. 



between p e and pi. What science can be gleaned from 
those excess events? 

We begin by providing a Bayesian perspective on this 
question. First, setting any threshold and discarding 
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events below that threshold amounts to "throwing away" 
data, and in the ideal case of arbitrary computing power 
and a perfectly known distribution of the detection statis- 
tic, discarding data weakens the analysis. So in the ideal 
case, p e would be set very low. In actual practice, how- 
ever, finite resources, an imperfectly known distribution 
of outliers, and diminishing infomation returns (for lower 
p e ) all will push data analysts to a value of p e not too 
far below p\. Probably the most important information 
encoded in the excess, sub-pi events is an improved esti- 
mate of the event rate 1Z (e.g., in units of Mpc _3 yr~ 1 ). 
As emphasized in Messenger and Veitch [60], a proper 
Bayesian analysis takes into account the value of p Cl and 
so always gives an unbiased estimate of 1Z. Including 
more events by lowering p e just "shrinks the error bars". 

Next we provide a crude, back-of-the-envelope argu- 
ment why, in the regime of rare detection, there will 
generally not be much information in the sub-threshold 
excess. In the regime of rare detection, the unique loud- 
est event will have SNR ~ 10. There are eight times 
as many events twice as far away, so we would expect 
of order 8 true events with 5 < SNR < 10. However 
the number of background events with SNR between 5 
and 10 will be of order 10~ 4 (the FAP at SNR - 10) 
times e -25 / 2 /e~ 100 / 2 , or ~ 10 12 (assuming a Gaussian 
distribution of noise events). The standard deviation in 
this expected number of events is 10 6 , which swamps the 
true excess. Of course, this calculation was for one (non- 
optimized) value of p e , but we believe that the general 
conclusion is robust: as p c is decreased below pi, the 
background rate increases far faster than the rate of true 
events, so we expect relatively little extra information 
from the sub-threshold events. 

In the next subsection we consider two non-Bayesian 
statistical excess tests currently used in LIGO searches, 
and show through simulations that, indeed, the sub-p! 
excess typically contains rather little useful information. 
We believe that a proper Bayesian analysis the sub-pi 
events would lead to the same qualitative conclusion, but 
we leave that calculation for future work. 



A. Results from Simulations 

In past LIGO searches, two statistical excess tests 
were employed: the Mann- Whitney- Wilcoxon U test 
and later the binomial test [5, 20]. However we are not 
aware of any systematic comparison of their detection 
power, so we compare them here. 

We ran a Monte-Carlo simulation in which each trial 
drew 20 GRB on-source candidates and, for each one, 10 6 
off-source trial loudest candidates. Off-source candidates 
were drawn from the S6/VSR2,3 distribution (see Sec- 
tion III B) . Each on-source candidate SNR was the max- 
imum of an off-source trial and a draw from a foreground 
p(p) oc p~ 4 distribution whose rate scaling could be var- 
ied. For each trial, we performed individual direct detec- 
tion searches on all 20 GRBs before performing a statisti- 




5 10 15 

Average number of direct detections 
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FIG. 2. Fraction of simulations in which the U and binomial 
tests detected a statistical excess in sub-threshold candidates 
vs the rate of direct detections. 



cal excess search on the set of remaining non-detections. 
Both direct detection of individual GRBs and statisti- 
cal excess detection had thresholds set by FAP = 10~ 4 . 
Figure 2 shows the fraction of trials in which the U and 
binomial tests were able to detect an excess versus the 
average number of direct detections among the 20 simu- 
lated GRBs. While the U test used all on-source and all 
off-source candidates, the binomial considered only the 
loudest four on-source candidates. 

We see our expectations confirmed: for both tests, the 
probability of detecting an excess in the sub-threshold 
populaton is never larger than several percent. Of the 
two tests, the binomial one is the more powerful. We see 
both tests become more powerful as the number of direct 
detections rises, until it reaches ~ 13. The decline in 
efficicency for > 13 detections is an artifact of the "rules" 
of our simulation in that there are always 20 true events, 
so as more of them exceed p\, fewer true events are left 
among the sub-threshold candidates. 



VI. SUMMARY AND DISCUSSION 

Based on LIGO-Virgo S6/VSR2,3 data and search pa- 
rameters for CBC-short GRB searches, we have quanti- 
fied how much more deeply into the noise one can dig 
with knowledge of the external trigger time; we find a 
coherent SNR threshold of 9.0 versus an all-sky coinci- 
dent threshold of 11.3 for a detection FAP of 10" 4 . [6, 7] 
have also estimated the reduction in the SNR threshold 
assuming Gaussian noise and making different choices in 
how to fold in EM information. The bleakness of [7] com- 
pared to our study is easily understood: if the measured 
distribution of SNRs falls off less steeply than a Gaus- 
sian, then reducing the search volume has a relatively 
larger effect on the detection threshold (at fixed FAP), 
and hence on the detection rate. While the Advanced de- 
tectors will not have the same distribution of high-SNR 
candidates SNRs as LIGO and Virgo did in S6/VSR2,3, 
the true distribution seems unlikely to be Gaussian, and 
the difference between our results and [ , ] provides some 
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measure of the importance of non-Gaussian backgrounds 
in assessing the value of triggered searches. 

Folding the threshold reduction into a large Mon- 
te-Carlo simulation, including the effects of short GRB 
collimation, general spin precession, and Advanced GW 
and GRB detector networks, we have estimated the rate 
of CBC-short GRB coincident detections. Assuming 
that all NS-NS systems produce short GRBs with a jet 
angle of 20°, we find that, relative to just an all-sky 
search, adding a search triggered by a Swift-like satel- 
lite increases the total number of CBC detections by 2%, 
but more importantly doubles the number of GW-GRB 
associations. A mission such as Fermi that has ~ 6 times 
the instantaneous sky coverage of Swift would contribute 
not quite an increase of ~ 12% to the total number of 
CBC detections, because the relatively poor sky localiza- 
tion would permit less of a reduction in the GW detection 
threshold. Although the calculated enhancement of de- 
tection rate is dependent on this and other assumptions, 
we believe it justifies the effort that is being spent on 
such triggered searches, given the extra scientific value of 
multi-messenger detections. 

The externally triggered GRB searches to date have 
attempted to detect a population of sub-threshold GWs. 
We performed simulations that show that typically there 
will not be a detectable excess until the rate of direct 
detections of individual sources is already high; hence 
it is highly unlikely that an excess-population test will 
provide the first strong evidence for CBCs. 
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Appendix A: Prior distribution of GW detector 
network's response to a CBC at a given distance 
and sky location 

In this section we discuss the distribution of a net- 
work's response to GWs from a source at a particular, 
known sky location and distance, whose orientation is 



unknown but whose inclination is restricted to be less 
than a maximum value t max = 8j . When studying an in- 
dividual GRB, we could treat this as a prior distribution 
for the strength of the signal received by a particular de- 
tector network. In the event of a non-detection, it would 
allow us to parameterize the excluded distance by the jet 
opening angle. 

Let the frequency domain GW strain received by de- 
tector X be h x (f) and the noise power spectral density 
of detector X be S x (/) . Defining the inner product 

( ^-^,^ ) -«B. J r 

and then the sum over all of the detectors, 



X\2 



x 



the coherent detection statistic p C oh 2 defined by Harry 
and Fairhurst [ 1 and the incoherent, coincident detec- 



tion statistic p c 



are noncentrally chi squared dis- 



tributed with the noncentrality parameter given by r 2 . 
(The coherent detection statistic has 4 degrees of free- 
dom whereas the coincident statistic has 2 degrees of 
freedom times the number of detectors.) We will first 
derive summary statistics of r: its minimum, maximum, 
mode, mean, and root-mean-square (rms). Then, we will 
derive the full distribution of r and study its qualitative 
features. 

Harry and Fairhurst [ ] introduce D, the luminosity 
distance of the source; Do, an arbitrary fiducial distance; 
the three angles (t, ip, <f>o) describing the orientation of the 
source, being respectively the inclination of the orbital 
plane to the line of sight, the polarization angle, and the 
orbital phase at coalescence; and the antenna factors of 
each detector, F? and F x , which are functions of sky 
location. They also define two waveform quadratures, ho 
and fin /2i which are nearly orthogonal such that 



and 



M/)l 2 
s x (f) 



4 Re 



df 



\K/2(m 

s x (f) 



X\2 



d/=(<r A ) 



K{f)K /2 {f) 
s x (f) 



d/wO. 



They define a further three quantities that combine the 
antenna factors of a network of detectors, 



x 

B = 5>*F*) 2 
x 

C = J2(a x F x )(a x F x ). 



X 



(Al) 



It is easily — though laboriously — shown that the detector 
response or noncentrality parameter r 2 depends only on 
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the antenna factors, distance, t, and ip, through 



1 A) 
8 D 2 



(A + B)(x 4 + 6x 2 + 1) 



+ ^(A-B) 2 +4C 2 (1 - x 2 ) 2 



cos a 



, (A2) 



where x = cosi and a = 4-0 — atan2(2C, A — B). 
For the purpose of concise parameterization of r 2 , we 
will also introduce xq — cost max , J 2 = A + B, and 
X 2 = y/(A — B) 2 + AC 2 /{A + B). J 2 describes the to- 
tal sensitivity of the detector network as the weighted 
sum-of-squares of all of the antenna factors. To lend in- 
terpretation to K 2 , we write Equation (Al) as 



H 



a^l a 2 F 2 a 3 Fl 
a^Fl a 2 F 2 <j 3 f} 



A C 
C B 



g 2 FI a 2 F 2 



a F 



Then, K 2 can be shown to equal the ratio of the difference 
of the eigenvalues of H to their sum. K 2 measures the 
extent to which the network is preferentially sensitive to 
just one polarization, in a way that is independent with 
respect to rotations of the detector network's coordinate 
system. If the network is equally sensitive to two orthog- 
onal polarizations, then K 2 vanishes. If the network is 
sensitive to only one polarization, then K 2 is unity. 
Last, we define u 2 as the distance-independent part of 



D 2 

W 2 [(» 



%x l 



1) +K 2 {1 - .T 2 ) 2 cosa] 



Our goal is to study the conditional probability density 
function (pdf), p(u | 0, (j), i < L max ), of u, assuming a fixed 



sky location and a censored orientation distribution. Un- 
der the assumption that the GRB emission is collimated 
within a an angle i max = 6j , the burst can only be seen if 
earth is placed inside this cone. As the probability is the 
same for random placement anywhere on the surface of 
this cone, the prior distribution of the inclination angle 
l between the line-of-sight to earth and the axis of the 
outflow is given by 



p(0 



sm l 



1 



where < i < 



(A3) 



since an area element on the cone is given by dA = 
sin i di d<f>. This represents a prior distribution on the 
direction of the system's orbital axis that is uniform in 
solid angle, restricted to polar angles < t max - 

Having parameterized the detector network's response, 
we proceed to derive the distribution of r for a given sky 
location and luminosity distance, as well as the mini- 
mum, maximum, mean, rms, and mode of the distribu- 
tion. Finally, as an example, we will apply our results to 
GRB 051103. 



1. Distribution and summary statistics 

a. Minimum and maximum 

The minimum value of u is obtained when the source 
is at the maximum inclination (t = t max or x = xq) and 
when cos a = — 1: 



W + 6x 2 + 1 - i^ 2 (l - xq 2 ) 2 . (A4) 



The maximum value is obtained when the source is at 
the minimum inclination (t ~ or x = 1): 



= J. 



b. Mean 



(A5) 



The mean response is given by 



x Jo 



J ^/(x 4 + 6x 2 + 1) +K 2 (1 - x 2 ) 
V8 2tt(1 - x ) 



cos a 



da Ax. 



It is possible to cast the integral over a into the form of a complete elliptic integral of the second kind, 

) da; (A6) 



/ ^/x i + 6x 2 + l + K 2 (l-x 2 ) 2 E I J . t ., , , , 

V2tt(1 - x ) J X( , \y x i + 6x 2 + l + K 2 (l-x 



2K 2 (l-x 2 ) 



where c. Root-mean-square 

E(k) = / \A - k 2 sin 2 a da. The rms response is 

Jo 
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1 ,2* j2 



x JO 



(x 4 + 6x 2 + 1) + K 2 (l - x 2 ) 2 cos a 
2tt(1 - x ) 



da dx. 



The integral over a kills the cos a term, leaving 
J 2 f 1 fx 4 + 6x 2 + l" 



u ln 



The integral over x gives 



1 - X 



dx. 



«ri 



(j27w) ^ x ° 4 + x ° 3 + llx ° 2 + 11x0 + 16 - 



(A7) 



and y = (x,it). By forming the Jacobian determinant 
\df/dx\, we find 



da; da 



1 

du du 
dx da 



dx du 



du 



da 



dx du. 



Using the inverse rule for derivatives, dy/dx 
(dx I 'dy)~ x , we solve for a as a function of u. 



a = arccos 



8u 2 /J 2 - (x 4 + 6x 2 + 1) 



d. Full probability distribution 
We have the prior distribution of (t, a), 

p(l, a) = - r where a G [0, 2?r] , l G [0, t max ], 

2tt(1 - cost max ) 

or equivalently, 

p(x, a) = — — r where a £ [0, 2tt], x G [xq, 1] . 

27r(l-a;o) 

To compute the conditional pdf of u, p(u | 9, </), i < t max ), 
the first step is to effect a change of variables from dx da 
to dx du. To do this, we write y = f (x) where x = (x, a) 



K 2 (l-x 2 ) 2 
and then differentiate with respect to u, 



da I6u 

du = ~ J 2 K 2 (1 -x 2 ) 2 



1 



8u 2 /J 2 - (x 4 + 6x 2 + 1) 
K 2 (l -x 2 ) 2 

16u 



[J 2 K 2 (l - x 2 ) 2 } 2 - [8u 2 - J 2 (x 4 + 6x 2 + l)] 2 



Now we may express the conditional pdf that we seek as 



p(u \Q,(j>,L< L n 







da 


r 


2 


du 



p(x, a) dx. 



The factor of 2 accounts for the two distinct values of a 
that give the same value of u. Altogether, 



p(u \9,(f>,l< L n 



16 



u/J 



KJ{1 Xq) J X1 j ^ R2{1 _ _ [g( u / j)2 _ ( x 4 + Q x 2 + ^ 



: dx 



(A8) 



We have to be a little careful about the limits of integra- 
tion; the quantity in the radical must remain positive. It 
has a zero at 



-3 + K 2 + ^8[l-K 2 +u 2 (l + K 2 )/J 2 } 



l + K 2 



The lower limit should be 



r 



The upper limit should be 



X2 



-3~K 2 +J8[l + K 2 +u 2 (l-K 2 )/J 2 ] ., , 
" T^K2 II IS. f= 1 

u/J iiK = l 

. (A10) 

Equation (A8) may be evaluated numerically using, for 
example, Simpson's rule. 



\ 



-3 + K 2 



+ ^8[l-K 2 + ^j 2 (l + K 2 )] 



e. Mode 



l + K 2 



The mode of the distribution occurs when the lower 
(A9) limit of integration ceases to clip against the minimum 
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value, where Neither w mC an nor the full pdf p(u \ 9,<f>,i< t max ) simplify 

much for the unrestricted inclination case. 



x = 
This occurs at 



-3 + K 2 + y/8[j - K 2 + m 2 (1 + if 2 )/J 2 ] 
l + K 2 ' 



^modc 



(-^j VV + 6V + i + if 2 (i-V) 2 . 

(All) 



2. Special case: unrestricted inclination 



If the inclination is unrestricted. t niax = 7r/2 or .r () = 0. 
then the above results simplify to 



Umode = J\/l + K 2 /V8 

Umin = J\/l-K 2 /VS 
^max J 
^rms 



= V2J/V5. 



3. Special case: one detector 



When the detector network consists of only one de- 
tector, C 2 = AB, K = 1, and J 2 = a 2 F 2 , where 



F 2 = F, 



. With these substitutions, 



"mode = ./(I + x 2 )/2 



"'min 
t*max 

Ltrms 



P(U \6,(j),L < L n 



Jx 

J 



(l + a WL^L| d, 



1G 



u/J 



7TJ(1 Xo) 7^ y p _ _ [g u 2/ J2 _ (j.4 + 6;E 2 + ^ 

X\ = y max (2u/J — 1, £o 2 ) 

X2 = ti/J. 



: da; 



4. Case study: GRB 051103 



GRB 051103 was an exceptionally short, hard, and 
bright burst detected by BETE, Suzaku, and Swift, and 
localized by IPN to an area consistent with the outer 
disc of M81 [23]. Owing to its brightness and hardness, 
a giant flare from an extragalactic soft 7 repeater (SGR) 
was a plausible progenitor. The Hanford 2-km detec- 
tor (H2) and Livingston 4-km detector (LI) were operat- 
ing at the time, so a targeted search of the GW data was 
undertaken. No candidate was detected, but the non- 
detection excluded a CBC event in M81 as the progen- 
itor [19]. Under the assumption that a CBC progenitor 
would have produced a collimated jet along the axis of 



strongest gravitational wave emission, Abadie et al. [19] 
placed 90% confidence lower limits on the distance of a 
CBC progenitor as function of jet angle. A collimated 
GRB in M81 was firmly excluded. 

As an example, we apply our distribution of detec- 
tor response to the problem of estimating the exclu- 
sion distance as a function of jet opening angle for 
GRB 051103. If we knew the GW search's detection effi- 
ciency for strictly face-on sources, then using our distri- 
bution for (u/J) we could directly calculate the excluded 
distance for any jet opening angle and any confidence 
level. Abadie et al. [ ] did not publish that detection ef- 
ficiency, but we can do a qualitatively similar calculation 
by extrapolating from their 90% exclusion distance for 
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FIG. 3. Prior distribution of (u/J), the distance-independent part of the detector response, normalized by the detector's 
root-sum-squared antenna pattern at a given sky location. The top row shows the probability density function itself, for jet 
opening angles of 15° to 90°. From left to right, the value of K 2 is varied from 0.25 to 0.1, smoothly varying from a detector 
configuration that has similar sensitivity to two polarizations to a configuration that is sensitive to only one polarization. The 
bottom row of plots shows five summary statistics of the distribution, the minimum, maximum, mean, mode, and rms. 



8j = i max = 10°, attempting to reproduce their exclusion 
distance at other jet opening angles. 

GRB 051103 occurred at 3 November 2011 09:25:42 
UTC. The H2 and LI horizon distances at this time 
for both a 1.4-1.4 M NS-NS event and a 1.4-10 M© 
NS-BH event are given in Table III, along with the an- 
tenna factors at this time in the direction of M81. For a 
NS-NS signal, this network has K 2 = 0.9601, and for a 
NS-BH signal, K 2 = 0.9602. H2 and LI had almost the 
same sensitivity up to a frequency independent factor of 
«2, so it is not surprising that the value of K 2 is almost 
the same for both the NS-NS and NS-BH signal models. 

In Fig. 4, we plot the 90% exclusion distance as a 
function of 0j ot from Fig. 3 of Abadie et al. [ ]. We 
have superimposed the mode of the detector response 
distribution, equation (All), scaled to match the pub- 
lished exclusion distance at 0j = 10°, as a dashed line. 
The value of (u/J) at which the cumulative distribution 
function (CDF) is equal to (1 — 0.9) is shown as a solid 
line. The inverse CDF agrees well with the NS-NS exclu- 
sion distance, but the mode agrees much better with the 
NS-BH exclusion distance than the inverse CDF. Ex- 
act agreement is not expected with either: as we have 
pointed out, a proper application to calculating exclu- 
sion distances would require knowledge of both the prior 
distribution of (u/J) and the sensitivity of the GW search 
to face-on sources as a function of signal amplitude. Fur- 
thermore, the analysis of Abadie et al. [19] includes a 



TABLE III. Antenna factors and horizon distances for H2 and 
LI detectors at time of GRB 051103. 



detector 


F+ 


F x 


-D-Hnsns 


DHnsbh 


H2 


-0.152 


-0.706 


8.3 


17.1 


LI 


0.348 


0.550 


18.1 


37.4 



Monte-Carlo integration over a range of masses, whereas 
our analysis fixes canonical choices of the masses. 



Appendix B: Enhanced GW amplitude of 
GRB-triggered sources 

In this paper, we assume that the emitted 7 rays are 
beamed within an angle (9b of the normal to the orbital 
plane and thus the binary inclination 1 must be less than 
6b- Where 6>b is small, we approximate the GW ampli- 
tude for GRB-triggered sources to be on average ~ 1.51 
times the isotropic detection-averaged amplitude for all 
CBCs, which is the instantaneous 1 = value. 

Using intermediate results from Appendix A, the az- 
imuthally averaged detector response to a binary whose 
orbital plane inclined at angle t relative to the observer's 
line of sight is proportional to 

(j) 2 ^(* 4 + 6z 2 + l) , (Bl) 
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FIG. 4. Exclusion distance as a function of i m ax = 0j- The 
filled triangles are the 90%-confidence exclusion distances 
from Fig. 3 of Abadie et al. [ I ! i] . The up-arrows (green in 
the online version) represent the NS-NS signal model, and 
the down-arrows (blue in the online version) the NS-BH sig- 
nal model. The solid curves show the 10% value of the in- 
verse CDF of (u/J), scaled to match the exclusion distance 
at i, max = 10°. The dashed curves show the mode of the dis- 
tribution of (u/J), also scaled to match the exclusion distance 
at t max = 10°. The pink band marks the distance to M81, 
3.63 ±0.14 Mpc. 



where x — cos t and J 2 = Fi + F 2 . The distance to 
which a GW source is detectable scales as u/J, so the 
number of detectable sources scales as (u/J) 3 . Thus the 
detection-averaged amplitude of the all sources that we 
observe within half-angle 6*b is 



A(x ) = 




where x = cos#b an d p{x) — (1 — x ) • 

Using Eq. B2, one easily shows that as long as the beam 
half-angle is < 25° , then average amplitude enhancement 
relative to an unrestricted distribution is at most a ~ 5% 
overestimate. 
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